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Abstract 

' In cosmological context, classical scalar fields are important ingredients for inflation models, 

, many candidate models of dark energy, symmetry breaking and phase transition epochs, and their 

■ consequences such as baryo and lepto-genesis. We investigate the formation of these fields by 

studying the production of a light quantum scalar field during the decay of a heavy particle, 
i' For simplicity it is assumed to be a scalar too. We discuss the effects of the decay mode, the 

, thermodynamical state of the decaying field, boundary conditions, and related physical parameters 

on the production and evolution of a condensate. For a simplified version of this model we calculate 
the asymptotic behaviour of the condensate and conditions for its contribution to the dark energy 
(N . with an equation of state close to a cosmological constant. We also discus the role of the back- 

reaction from interactions with other fields and expansion of the Universe on the evolution of the 
condensate. 
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^ ■ 1 Introduction 
O 

• Symmetry breaking, phase transition, and related phenomena such as appearance of a dynamical mass 

I (Higgs mechanism), superconductivity, inflation, leptogenesis, and quintessence field as a candidate 

Q^' for dark energy in early universe and cosmology are based on the existence of a classical scalar field. 



As the physics of the Universe and its content in its most elementary level is quantic, this scalar field, 
i) ', fundamental or composite, is always related to a quantum scalar field. 

A classical field is more than just classical behaviour of a large number of scalar particles. In a quantum 
system, fields/particles are in superposition states and are quantum mechanically correlated with each 
I others. Decoherence process removes the superposition/correlation between particles. However, this 
' does not mean that after decoherence of scalar particles, they will behave collectively like a classical 
field. A simple example is the following: 

Consider a closed system consisting of a macroscopic amount of unstable massive scalar particles that 
decay to a pair of light scalar particles with a global SU{2) symmetry and negligible interaction. If 
the unstable particle is a singlet of this symmetry, the remnant particles are entangled by their SU (2) 
state. After a time much larger than the lifetime of the massive particle, the system consists of a 
relativistic gas of pair entangled particles. If a detector measures this SU{2) charge without significant 
modification of their kinetic energy, the entanglement of pairs will break i.e. the system decoheres 
and becomes a relativistic gas. The equation of state of a relativistic ideal gas is Wrei = P/p ~ 1/3- 
In contrast to an operator in quantum field theory, a classical scalar field ip{x) is a C-number. Its 
density py,, pressure and kinetic energy are defined as: 

p^^K^ + V{^) (1) 
P^^K^- V[^) (2) 

= ^g^'d^^d,^ (3) 

V{(p) is the potential presenting the self-interaction of the field '^{x). When it is much smaller than 
kinetic energy , one obtains P^ ^ p^. On the other hand, if V{ip) » K^, « —Pf Therefore 
in general, a relativistic gas and a scalar field don't share the same equation of state. Thus, the proof 



of decoherence in a system is not enough when a classical scalar field is needed to explain physical 
phenomena. 

According to canonical quantization procedure, classical observables are replaced by operators acting 
on a Hilbert or Fock space of states respectively for a single particle and for a multi-particle quantum 
system. The expectation value of these operators are the outcome of measurements. Therefore, it is 
natural to define the classical observable related to a quantum scalar field as its expectation value: 



where \ ^) is the state of the quantum system - an element of the Fock space. In analogy with particles 
in the ground state in quantum mechanics, the classical field ip{x) is also called a condensate. Using 
canonical representation it is easy to see that for a free quantum scalar field (^'|<1>|^') = 0. Therefore, 
a necessary condition for the buildup of a classical scalar field is an interaction (see also Appendix R|) . 

In quintessence models a classical field is the basic content of the model and the source of the dark 
energy. Although in the framework of popular particle physics models such as supersymmetry, su- 
pergravity, and string theory many works have been concentrated on finding candidate scalar field to 
play the role of quintessence [1], little effort has been devoted to understand the necessary conditions 
for a quantum scalar fields to condense in a manner satisfying very special characteristics needed for a 
quintessence field. For instance, such a condensate must have a very small density, much smaller than 
other content of the Universe (smallness problem) . Present observations show that dark energy has a 
behaviour close to a cosmological constant i.e. with expansion of the Universe either its energy density 
does not change or varies very slowly. Such a behaviour is not trivial. In the classical quintessence 
models usually the potential or other characteristics of the models are designed such that a tracking 
solution is obtained. However, apriori it is not trivial to prove that in the early Universe a quantum 
field could produce a condensate with properties similar to dark energy. 

The purpose of the present work is to fill the gap between quantum processes producing various 
species of particles /fields in the early Universe - presumably during and after reheating - their classical 
component as defined in (j4]). In another word we want to see how the microscopic properties of matter 
is related to macro-physics and vis versa. As the quantum physics of that epoch is not well known, 
we consider the simple case of a scalar field - quintessence field - in interaction with two other scalar 
fields as a prototype process, and study the evolution of the classical component (condensate) of the 
quintessence field. Between many possible types of quantum scalar field and interaction models, we 
specially concentrate on a class of models is which the scalar field is one of the remnants of the decay 
of a heavy particle. Motivations for such a model are the results of studying the effects of a decaying 
dark matter on the equation of states of the Universe [2]. It has been shown that a FLRW cosmology 
with a decaying dark matter and a cosmological constant at late times behaves similar to a cosmology 
with a stable dark matter and a dark energy component with w = P/p ^ — 1. This is effectively what 
is concluded at least from some of present supernovae observations [3]. Recently, the same effect has 
been proved for the general case of interaction between dark matter and dark energy [4] . It has been 
also shown [5] that if a decaying dark matter has a small branching factor to a light scalar field, this can 
explain both observed density and equation of state of the dark energy without extreme fine-tuning 
of the potential or coupling constants. In other words, such a model solves both the smallness and 
the coincidence problems of the dark energy. These studies however are based on the assumption of a 
classical scalar field (a condensate). The present work should complete this investigation by studying 
the formation and evolution of classical component from quantum processes. 

In Sec. [2] we construct the Lagrangian of a decaying dark matter model. We consider two decay modes 
for the heavy particle and use the closed time path integral method to calculate the contribution of 
interactions to the condensate. The same methodology has been used for studying infiation models [6], 
late-time warm inflation [7], the effects of renormalization and initial conditions on the physics of 
infiation [8], baryogenesis [9], and coarse-grained formulation of decoherence [10]. In Sec. [3] we obtain 
an analytical expression for the asymptotic behaviour of the condensate and discuss the importance 




(4) 



of the back-reaction of the quantum state of the Universe on the evolution of the condensate and 
properties of the dark energy. We summarize the results in Sec. [H In Appendix [A] we make a remark 
about relation between decoherence and formation of a condensate. In Appendixes [B] we show that 
the effect of a non-vacuum state on the Green's function can be included in the boundary conditions. 
Appendix O presents the solution of the evolution equation of the field in matter dominated era. 

2 Decay in an Expanding Universe 

We consider a simple decay mode for a heavy particle X with only 2 types of particles/fields in the 
remnants: a light scalar $ - light with respect to the decaying particle - and another field A of an 
arbitrary type. In fact, in a realistic particle physics model, most probably A will not be a final stable 
state and decays/fragments to other particles. Therefore, it should be considered as an intermediate 
state or a collective notation for other fields. In the simplest case studied here all the particles are 
assumed to be scalar. Extension to the case where the decaying particle X and one of the remnants 
are spinors is straightforward. The quintessence field $ must be a scalar to condensate. Nonetheless, 
in the extreme densities of the Universe after reheating, apriori the formation of Cooper-pair like 
composite scalars is possible if the interaction between spinors is enough strong. Thus, $ can be such 
a field, but for the simple model studied here we ignore such complexities. 

The simplest decaying modes are the followings: 



Diagram ([5]-a) is a prototype decay mode when X and $ shares a conserved quantum number. For 
instance, one of the favorite candidates for X is a sneutrino decaying to a much lighter scalar field 
(e.g. another sneutrino) carrying the same leptonic number [13] [H]. With seesaw mechanism in the 
superpartner sector (or even without it [U]) if SUSY breaking scale is lower than seesaw scale, a mass 
split between right and left neutrinos and sneutrinos will occur. As the right-hand neutrino super-field 
is assumed to be a singlet of the GUT gauge symmetry, it has only Yukawa-type of interaction. In 
such a setup X can be a heavy right sneutrino decaying to a light sneutrino with the same leptonic 
number and a pair of Higgs or Higgsino [15] . In place of assuming two A particles in the final state we 
could consider them as being different A and A' . But this adds a bit to the complexity of the model 
and does not change its general behaviour. For this reason we simply consider the same field. Diagram 
([5]-b) is representative of a case where X and A are fermions, or $ carries a conserved charge |16j . 

It was easier to consider a simple 3- vertex X ^ ^ + A similar to what is considered in Ref. [?]• 
However, such a vertex does not always allow simultaneous conservation of energy. Moreover, as 
we will see later, more complex diagrams considered here will show how the decay mode affects the 
evolution of the condensate component. 
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The corresponding Lagrangians of these effective interactions are the followings: 



i:-x 



d^Xy 
d^Xy 

d^Xy 



2 ^ 2 ^ n' 



-ig^XA^, For ©-a 
^ yg^'^XA, For (P-b 



(6) 
(7) 
(8) 
(9) 



In addition to the interaction between X, $ and A we have assumed a power-law self-interaction for <I> 
and A. If A is a collective notation for other fields in the actual model, its self-interaction corresponds 
to the interaction between these unspecified fields. Again for the sake of simplicity in the rest of this 
work we consider A' = 0. The unstable particle X is assumed to have no self-interaction. 

Although the model presented here is quite general, for physical and observational reasons we con- 
centrate on the case of a heavy X particle as a candidate for the dark matter, $ as a quintessence 
field, and interactions in ([5]) as candidate interactions for the decay of the dark matter and production 
of what is observed as dark energy. It is therefore necessary that X and $ have only a very weak 
interaction. Therefore couplings A and g must be very small. 

In a realistic particle physics model, renormalization as well as non-perturbative effects can lead to 
complicated potentials for scalar fields. An example relevant for dark energy is a pseudo-Nambu- 
Goldston boson as <I> and potentials with a shift symmetry [T7]. These models are interesting for the 
fact that the mass of the quintessence field does not receive quantum corrections and can be very 
small. Moreover, they can be easily implemented in SUSY theories in relation with right-neutrinos 
and sneutrinos (as candidate for X). The power-law potential considered here can be interpreted as 
the dominant term in the polynomial expansion of the potential. In any case, general aspects of the 
analysis presented here do not depend on the details of the particle physics and self-interaction, and 
can be applied to any model. The solution of the field equation and numerical quantities are however 
sensitive to the particle physics. The purpose of the present work is to investigate the behaviour 
of this model and to find what is most important for the formation and evolution of a condensate 
with characteristics similar to the observed dark energy. We leave the application of this analysis to 
realistic particle physics models to a future work and consider only simplest cases which are analytically 
tractable and permit exact or approximate analytical solutions. 



We decompose ^{x) to a classical (condensate) and a quantum component: 

^{x) = ip{x) + (l){x) {^) = = ip{x) (0) = (^l^l^r) = 



(10) 



Note that in (jlOp both classical and quantum components depend on the spacetime x. In studying in- 
flation it is usually assumed that very fast expansion of the Universe washes out all the inhomogeneities 
and the condensed component is homogeneous. As we are studying the evolution after inflation, the 
distribution of unstable X can have non- negligible inhomogeneities, specially if the decay is slow and 
perturbations have time to grow. 

We assume {X) = and {A) = 0. Justification for these assumptions is the large mass and small 
coupling of X which should reduce their number and their quantum correlation. In other words, when 
mass is large, the minimum of the effective potential for the classical component is pushed to zero (see 
and (fT3]) below). We find a quantitative justification for negligible condensation of massive fields 
in Sec. Ql 



The Lagrangian of $ is decomposed to: 



/ 



d —g 



(gipXA^ + g^XA^ For ©-a 

|gV?2j^A + 2g^(l)XA + g(/>2XA For ©-b 



A / n 
n V ^ 



+ 



(11) 



Lagrangians and are the same as ^ with respectively and $ ^ After replacing 

quantum terms by their expectation values, this Lagrangian leads to the following evolution equation 
for the condensate component: 



^d^.i^g^-'d.y,) + m|<^ + ^ ^(i + 1) . ^ vHr-'-') - g{XA^) = 

n \i + IJ 

For ©-a (12 

, n—l 
A 



1 A / \ 

-^^^(y^^^'^S.v^) + + - ^(i + 1) . "I 



i iin—i—l\ 



2g^{XA) - 2g{(t)XA) = 

For ®-b (13) 



In the case of the interaction mode ©a), the interaction Lagrangian depends linearly on <I> and appears 
as an external source in the field equation (jl2|) . For both decay modes, if n ^ 2, the term i = in the 
sum of the self-interaction terms also contributes to the non-homogeneous component, and the term 
i = 1 contributes to the effective mass of the classical field (p. Note that what we call non-homogeneous 
component or external source terms have in fact implicit dependence on ip . The reason is the coupling 
between quantum interactions and the evolution of the classical component. We will show later that 
these terms play the role of a feedback between production and evolution of the condensate. In fact, © 
b) has a reach structure and various evolution histories are possible depending on the value and sign of 
g, the coupling to X, self-coupling A, and the order of the self-interaction potential n. For instance, the 
mass can become imaginary (tachyonic) even without self-interaction leading to symmetry breaking. 
Tachyonic scalar fields have been suggested as quintessence field specially in the framework of models 
with w < —1 [18]. The decay mode ©a) by contrast has a field equation very similar to the classical 
model studied in Ref. [5] . We discus in detail the differences of these decay modes in the next sections. 

In addition to interaction terms in the Lagrangian of the total field the Lagrangian of the purely 
quantic component (p includes terms depending on the classical component ip. Derivative term 
^g^'^ {dp^ipdu4) + d^4)di,ip), mass term TJi^ipcj), i = n — 1 term in the self-interaction, and 2gip(pXA 
are linear in (p and only affect the renormalization of the propagator [8] . For n ^ 2 the term i = n — 2 
in the self-interaction sum contributes in the effective mass of cj) and makes it time dependent. 

We use Schwinger closed time path (also called in-in) formalism to calculate expectation values. 
Recent reviews of this formalism are available [11] and here we only present the results. Zero-order 
(tree) diagrams for the expectation values (fT2l) and (fT3]) are shown in (fT4l) . (fTSi) and (fTBl) . The 
next relevant diagrams are of order g^ and negligible for the dark energy model. One example of 
higher order diagrams is shown in (I14p . These types of diagrams are specially important for studying 
renormalization in the context of a realistic particle physics model. Thus for the phenomenological 
models considered here, we can ignore them. 





(14) 



g^\XA) 




+ ... 



(15) 
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Diagrams for self-interaction terms in (jl2p and (jl3p are similar to (jl5p with n — 1 external lines of 
type (j) or ip. The dash line presents the classical component of <1>. The corresponding expectation 
values at zero order are: 



g^{XA^) 
g^''{XA) 
9^{<I^XA) 



gip{x 
gip{x 



-gd y 



G>{x,y)G>{x,y)G>{x,y)-G<{x,y)G<{x,y)G<{x,y) 



-gd^y 



G> (x, y)G> (x, y) - G< (x, y)G< (x, y) 



-gdS 



G> (x, y)G^ (x, y)G> {x,y) - G< (x, y)Gi (x, y)G< (x, y) 



Future and past propagators G^ and are defined as: 



G>{x,y) 
G<{x,y) 



-i{'tp{x)'tp\y)) = -itT{^p{x)^^y)p) 
^i{'tp\y)'>p{x)) =^itr{^p^y)'il){x)p) 



(17) 
(18) 
(19) 



(20) 
(21) 



where tp{x) presents one of (p, X or A fields and p = |^')(^| is the density (projection) operator for 
the state |^). The upper and lower signs in (j2ip are respectively for bosons and fermions. Definitions 
(|20p and ()2ip correspond to the general case of a complex field. Here we only consider real fields and 
therefore ■ip{x) = -(/'^(x). Feynman propagators are related to G^{x,y) and G^{x,y): 



GF{x,y) ^ -i(rV^(x)V^t(y)) = G'>(x,y)e(xO-yO)+G<(x,y)e(yO-xO) 
GF{x,y) ^ -i(fV(x)VHy)) =G'>(x,y)e(y°-xO)+G<(x,y)G(xO-2/0) 



(22) 
(23) 



The next step is the calculation of propagators. 



2.1 Propagators in an expanding universe 



Feynman propagators Gp{x,y), i = (p, X, A can be determined using field equations from La- 
grangians ([6])- ([9]). As the concept of Green's functions is only applicable to the linear differential 
equations, we have to linearize the field equations and treat interactions perturbatively, i.e. as quan- 
tum corrections. Free propagators of (j), X and A are: 



--d^iV^g'^'d^GUx - y)) + (mi + (n - l)A<^"-^)G^(x - y) 



-M^gg^"d,G'F{x - y)) + mfG^pix - y) 



.(5^(x - y) 



.5^{x-y) . 
-I ^= — I 



(24) 
X,4;25) 



The free propagator of (j) is independent of the type of interaction with X and A and therefore, 
equation (p^ is valid for both interaction models presented in ([5]). Note also that G^(x — y) is 
coupled to the classical field (p even at the lowest quantum perturbation order. On the other hand, 



evolution equations (|T2|) and p3]) . depend on the interaction between quantum fields cj), X and A. 
Therefore, this model is coupled at all orders. 

To proceed, we neglect the effect of spatial anisotropy and consider a flat homogeneous metric with 
synchronous or conformal time: 

ds^ = dt^ - a^{t)5ijdx'dx^ = a^{'n){drf - dijdx'dx^) dt = adr] (26) 

where t and rj are respectively comoving and conformal time. It is more convenient to write evolution 
and propagator equations with respect to conformal time rj. After a variable change: 

X = aip T = a(j) l = aX k = aA (27) 

and by using the metric (|26p . the evolution equation of the classical field x takes the following form: 



For dSD-a (28) 



X" - S^A^JX + {a'rnl - -)x + + 1) • ! i xHT^"'^') - 2gx(XA) - 2g(TXA) = 

For (P-b (29) 



i=0 



Propagator of quantum fields T, X and A are: 



^Gj{x,y) - 6,,did,Gj{x,y) + {a'ml -- + (n - l)Aa V""')GJ(x, y) = - /'^^ ^^ (30) 

I^G^x, y) - S^jdidjG'^ix, y) + (a^m^ - ^)G%{x, y) = -i^-l^^-A i = x, A (31) 
Gl = a{ii)G^p , G^ = a(r/)G^ , = a(r/)G^ (32) 



with /' = df /dr]. Note that the classical component of $ appears as a spacetime dependent mass term 
for its quantum component (p (or equivalently T). Fourier transform can be applied to (|3ip . but the 
spacetime dependence of coefficients in (j28p . (|29p and (j30p makes this method useless for solving these 
equations. Nonetheless, if these terms are small and/or vary slowly, one can first ignore them and 
solve the equations. Then, by applying the WKB method, a more precise solution can be obtained. 
Equations (j30p - (j3ip are second order partial differential equations, and therefore propagators are linear 
combination of two independent solutions of the associated homogeneous equations. The delta function 
on the right hand side however leads to a discontinuity which appears as a consistency condition for 
the solutions and fixes the ambiguities in the propagator solution. This will be described in details in 
the rest of this section. 

X particles are presumably produced during reheating epoch |12j and begin their decay afterward. 
In this epoch relativistic particles dominate the density of the Universe. Thus, we first consider 
this epoch. Fortunately, for this epoch the homogeneous field equation has an exact solution. In 
matter dominated epoch only for special cases an analytical solution exists. They are discussed in the 
Appendix [Cl The expansion in the radiation domination epoch has the following time dependence: 

a = ao(^] =ao— (33) 



And thus a" = 0. After taking the Fourier transform of the spatial coordinates and neglecting the 
(/3-dependent term, the solutions of the associated homogeneous equation of ([30]) and (j3T|) are well 



known [l9] [20]: 
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d^xU\x)e 



ik.x 



$ , X , A 



Oii 



:i+i)\/B,, 



1 + 



m, 



2 - alHoTTli 



(34) 

(35) 
(36) 

(37) 



where qq and -ffo are respectively expansion factor and Hubble constant at the initial conformal time 
770 • The function Dq{z) is the parabolic cylindrical function. From now on for simplicity we drop the 
species index i unless when its presence is necessary. We call two independent solutions of in a 
general basis and Vfc. If we want these solutions correspond to the coefficients of the canonical 
decomposition of (p ( (j66p in Appendix [B]) . we must choose a basis such that Vk = ^l- ^'^^^ 
this work we only consider this basis. The corresponding equation for the free Feynman propagator 
(free 2-point Green's function) is: 



6{ri — rj' 



(38) 



When rj ^ rj', (j38p is the same as (j34p and therefore solutions of the former is a linear combination of 
two independent solutions of the latter. According to the definition of Feynman propagators (j20p and 
([2T]) it can be divided to past and future propagating components and G^. With 77 <-> ry' these 
propagators change their role: ^ . Therefore, G{rj,rj') has the following expansion: 



iG{r],r]') 



A<Uk mi (V) + [B<Ul iv)Uk iv') 



where , , A'j: and B"^ are integration constants. In the Appendix! 



(39) 

I we show that for the free 
propagators i.e. at the lowest perturbation order, if the state |^') is not vacuum, it is possible to 
include its effect in the boundary conditions imposed on the propagator. Comparing ([39]) with ()70p 
in the Appendix [Bj the relation between these constants and the initial state can be concluded: 



Skk,\^kik2...kj'^ 

It is easy to see that with these relations the consistency condition defined as: 



E E 

i kik2...kn 



G>{v,v'] 



G<{i],i]'] 



r)=r)' 



(40) 
(41) 



(42) 



is automatically satisfied. Therefore propagators over a non-vacuum state ^ only depends on this state 
and the solutions of the field equation considered as the free particle states. On the other hand, these 
solutions depend on two arbitrary constants and d\. which should be fixed by initial conditions. We 
have used decomposition (|39p along with canonical decomposition of the quantum field (see (j66p to 
()69p in Appendix B)to fix the integration constants and d^. The advantage of this method is that it 
explicitly relate dynamical constants to the physical properties of environment in which the quantum 
field is living. 

There is one more consistency condition that propagators should satisfy. By integrating two sides of 
equation (j38p with respect to rj in an infinitesimally region around r]' we find the following constraint: 



u^mm-Ukmk*{v) 



a{r]) 



(43) 



This relation fixes one of the dynamical constants in ()36|) . It can be chosen to be the normalization 
of the propagator. It is however interesting to note that multiplying both sides of (j43p with an 
arbitrary constant rescale a(ry) which is equivalent to redefinition of oq. Rescaling of oq is equivalent 
to redefinition of coordinates and therefore is not an observable. In Minkovsky spacetime the scale 
factor a is fixed to 1, and therefore there is no place for rescaling. Thus in a Minkovsky space the 
normalization of the propagators is an observable and affecting final results. The scaling properties in 
FLRW or De-Sitter metric is a consequence of diffeomorphism invariance in the framework of curved 
spacetimes and general relativity. 



2.2 Initial Conditions 

Field equations are second order differential equations and need the initial value of the field and 
its derivative or a combination of them to be totally described. The general initial conditions for a 
bounded system, including both Neumann and Dirichlet conditions as special cases, are the follow- 
ings [21]: 

ni'df.U = iKU , ni'd^JA* = -ilCU* , gf^yii^rf = 1 (44) 

The 4-vector is the normal to the boundary surface. If the boundary is space-like, the normal n'^ 
can be normalized to o~^(l, 0, 0, 0), then: 

a-^d^U = iKU , a'^d^U* = -ilC*U* (45) 

Constants fCi and /C/ depends on k. In a general boundary problem the boundary conditions must 
be defined for all the boundaries. Thus, in a cosmological setup the initial conditions (j45p must be 
applied to a past (initial) and future (final) boundary surfaces. This is the strategy suggested in 
Ref. |21j . The past and future boundary conditions are respectively applicable only to past and future 
propagators. Assuming different values for K, on the past and future boundary, one finds: 



i4M 

' ajUk{r]j) 



Indexes i and / refer to the value of quantities at initial and final boundary conditions. These boundary 
conditions relate /Cj and /C/ to Ck and dk in p6|) . In a cosmological context although fCf apriori can 
be decided based on observations, the value of /Cj is unknown and leaves one arbitrary constant in 
the solution or it should be selected according to a special model considered for the physics of early 
universe. This arbitrariness of the general solution, or in other words the vacuum of the theory, is well 
known [23]. In the case of inflation, this leads to a class of possible vacuum solutions called a- vacuum. 
For instance if: 



/C, = /C/ = ^A:Va2 + m2 (47) 
one obtains the well known Bunch-Davies solution 1211. 



Another way of proceeding is using the consistency condition (j43p to fix one of the arbitrary constants 
and applying the boundary condition (|46p only to one of the initial or final 3-surfaces. Although 
this does not solve the problem of arbitrariness of IC and its k dependence, it reduces it to only one 
of the boundary surfaces, for instance to the final 3-surface, and make the choice of (I47p physically 
motivated. Besides, in fixing only one of the boundary conditions, the causality of the solution is 
transparent - the state of the second boundary is directly related on the choice of the first one through 
the evolution equation. 

Another and somehow hidden arbitrariness in this formalism is the fact that apriori k dependence 
of the boundary constant IC does not need to be the same for all the fields of the model. However, 
different k dependence for the boundary conditions breaks the Equivalence Principal. Similarly, a 
value different from (j47p for IC will lead to the breaking of the translation symmetry [21] [8] . In the 
context of quantum gravity the violence of both of these laws are expected and therefore, in a general 
framework they should be considered. 



2.3 WKB approximation and back reactions 



Finally after finding the solution of linearized field equation ([3l|) and corresponding propagator ([38|) 
for the field T, we should add the effect of spacetime dependence of the mass term. According to the 
WKB prescription we must replace a^rj in U with: 



where (fk is the Fourier transform of ip{x) and q$ is defined in (j37p . With this correction apparently 
we have the solution for all the propagators at lowest order. However, <p{x) (or equivalently x) evolves 
according to equation ([28|) or ([29|) that depend on the expectation values of the quantum fields, 
specially (p. Therefore the propagators and the condensate are coupled even at lowest order, and only 
through a numerical calculation a final solution for each of them can be obtained. In fact the reason 
for coupling of classical and quantum fields is the self-coupling of $. If the self-interaction of $ is 
negligible, the solutions of evolution equation of the classical component (I28p or (I29p are similar to 
the field equation of the quantum component. In presence of a self-interaction however this equation 
is non-linear and must be solved numerically. 

In addition to the self-coupling, the state I'l') for which we calculate propagators and expectation 
values is a source of back-reaction (classical effects of back-reaction are studied in Ref. [22] ) . It defines 
the quantum state of the Universe at the time when the interactions or more exactly the decay of X is 
studied. However, |^) evolves due to interactions between species and the expansion of the Universe. 
They are responsible for the variation of the number of particles and their momentum distribution 
and therefore evolution of I'l'). In fact, in the path-integral formulation of closed path integral, a 
non- vacuum state |^) adds a functional integral to the path-integral [TTj that presents the projection 
of the state on a predefined basis. When the Fock space is evolving, as it is the case in the cosmological 
context, the projection of |^) evolves and path integrals become unfactorizable. 

As an example we consider an initial non-zero distribution for the X particles and no ^ or $ particles. 
At a later time, the number of X particles in |^') is reduced due to their decay and the expansion of 
the Universe. Their temperature or mean kinetic energy if they don't have a thermal distribution also 
decreases. On the other hand, a non-zero number of <I> and A particles are created. The latter at their 
production are relativistic and non-thermal. But if they have self- interaction and/or interaction with 
other fields which we ignored here, their momentum distribution will change both by interactions and 
due to the expansion of the Universe. This means that l^*) will change and its variation is refiected on 
the evolution of the classical component (p, the quantum component (p, the expansion factor a{r]), and 
the thermalization of A which we assume to have interactions with other particles (see also Appendix 



3 Evolution of the classical field without self-interaction 

To get an insight into the evolution of the classical component if, in this section we neglect self- 
interaction of <I> which is the main source of the non-linearity and coupling of the equations. We find 
an analytical solution for the evolution of condensate and discuss the difference between two decay 
modes, as well as the effect of the other parameters. The evolution of |\I') is introduced by a simple 
parametrization. We discuss its effect and determine the range in which the condensate can have a 
behaviour similar to the dark energy. 




(48) 



ED. 



3.1 Expectation values 



Neglecting the self- interaction of expression ([36|) is the exact solution of the field equation and ([39] 
is the exact propagator. Therefore, we can calculate the expectation values p!7|) to (fT9]) : 



g{XA){x) 



g{TXA){x) 



-^9 



Gl^iv,v)Gi>{v,v')Gl>iv,v') - Gifiv,v')Gi<{v,v')Glf{7i,v') 



Gi^{v,ri')Gl>{v,v') - GifM)Glfiv,v') 



d^kid^k2d^kze~'^-^^^+^^+''--'^5^^\ki + k2 + kz) / d-q' ^ 
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(49) 
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The most important aspect of these integrals for us is their time dependence because these expectation 
values contribute to the build-up and the time evolution of the condensate. Although the spatial 
spectrum is important specially for observational verification of the model, its role is secondary and 
comes later. 



Regarding 
following: 



and (|39|) . it is easy to conclude that time integrals in (|49|) - (j51|) are similar to the 



j=i 



N 



1 ± 



z = aiT]' (3i G {±, ±i} 7i 



i=l 



(52) 

where and Bi are defined in ()37p and N is the number of fields in the expectation value brackets. 
For simplicity we use the product index in (j52p for distinguishing the field species too. The constant 
coefficients of these integrals are of the form A{CiCl* with e {A^> , Bi> , Ai< , Bi<} , G {4,4} 
and C" G {4*'^*} corresponding species. These integrals don't have analytical solutions. 

Using the asymptotic properties of the parabolic cylindrical functions, we estimate the late time 
behaviour of these integral^: 



~ e * z^ 



^ _ g(g- 1) ^ 



Dq{z) ~ 2'ie- 
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(53) 
(54) 



The validity of these regimes depends on the mass of the corresponding field and on the cosmological 
parameters at the initial time r/o: 



\l3iZQ\ 



2nii 



(55) 



For a light its mass at the initial time tjq can be comparable or lighter than the Hubble constant 
Hq and therefore approximation ()54p is applicable. As for X and A, it seems unlikely that in any 
relevant particle physics model for these fields be small and therefore approximation (|53p must 
be applied. 



'^Validity of both approximations depends on the spatial scale ki and the constant Bi for each species. The initial 
time r]o and Hubble constant Hq also play important roles in these approximations and in the general behaviour of the 
model. Here we take these conditions for granted. 



Applying these approximations to ()52p . at lowest order in r] this integral has the following time 
dependence: 



where A^' is the number of fields to which approximation (|53p can be applied. Index i refers to all the 
fields. Indexes i' refers to fields for which approximation (153p is applicable and j for ones with (j54p 
approximation. If j = 0, the nominator in (j57p is 1. In this case the effect of the fields for which in 
the integration interval is small contributes only in the oscillating term of (j56p and in the higher-order 
terms of the polynomial expansion. 

For estimating the expectation values (I49p -()5ip in addition to time dependence included in integral 
(|52p one has to take into account the time variation of the state |^') both due to the expansion of the 
Universe and due to the back-reaction of the decay on the density of particles. Neglecting the effect 
of the decay on the density for a short duration after the massive production of X particles during 
reheating, the main source of the time evolution of j^') and \'^kik2...k„\'^ is the expansion. Assuming 
a very heavy X, it can be considered as non-relativistic at the time of its production and its density 
decreases by (x r]^. The other two fields A and specially $ are relativistic and their density decreases 
by a factor of oc r/^. 

As the lifetime of X particles is very long and thus g is very small, the contribution of A and (j) particles 
to the total number, energy density, and entropy of the relativistic matter at this epoch is very small. 
Therefore, one expects that after an initial fast increasing trend of the expectation values in (I49p -()5ip. 
they will slow down and probably approach a constant or begin to decrease. The turning point depends 
on the expansion rate, lifetime of X particles or equivalently the coupling constant g, the density of 
X particles at the end of reheating through its effect on |'I'fcifc2...fcnP' finally the mass of $ and 
A particles that can be estimated from (I56p . Larger the number of fields with m/Ho ^ 1, slower 
the time evolution of (j56p and sooner the turning point of the expectation values. This observation 
is consistent with our initial assumption of no condensation for X and A. We expect that only for 
the field $ there is a time interval in which m^/HQ < 1. Thus, at late times (XA^)(x) decreases faster 
than (XA)(x) and (TXA)(x). In addition to the time dependent mass term due to (XA^)(3;) in decay 
mode 5, the time evolution of expectation values in the two decay modes a and b are different, and 
therefore the condensation evolution for these modes are different too. Nonetheless there can be an 
exception to this argument. If we consider j4 as a single quantum field and not a collective notation, 
it is conceivable that at very early time symmetries keep it massless and it receives a dynamical mass 
after an epoch of symmetry breaking and phase transition. In this case the classical field x will have 
a faster gross and at later times, a sudden change in the time evolution of the expectation values 
()49p -(f5T]) and x is expected. Note that this argument is based on the back-reaction of the X particles 
decay on |^). 

In the next subsection we show that the most important term in the late time evolution of x is the 
non- homogeneous term in (j28p and ()29p . The coupling g appears as a time independent constant 
and changes the amplitude of x- This means that the density of the condensate is proportional to 
g^ (X T. This confirms the results of the classical treatment of this model in Ref. [5]. Regarding 
the coincidence and smallness problem of the dark energy, the close relations between mass, number 
density and lifetime of X and the mass and amplitude of the condensate is an evidence for an intrinsic 
feedback between dark matter and dark energy in this model. Moreover, this shows that the density 
of the dark energy is not more fine-tuned than the mass difference between left and right neutrinos 




(56) 



F{k) oc 



(57) 



in seesaw mechanism. In fact this similarity hints to a seesaw-hke mechanism for the large difference 
between the mass of X and ^ and supports the idea of right neutrino/sneutrino for X. 



Although here we are mostly interested in the decay of a meta-stable heavy particle, it is also interesting 
to see how a light scalar produced during the fast decay of another field, e.g. inflaton, can condensate 
and whether this condensate can last for long time and contribute to the build up of the dark energy. 
As before, in this setup the expectation values (|l9]) - ([5T]) will have an initial rise. But, depending on the 
lifetime of X particles is shorter or longer than the turning point time scale, either the latter happens 
earlier and x rapidly changes its slope from power-law to exponential decrease soon after the turning 
point, or power-low behaviour lasts for much longer time and the exponential break happens late in 
the lifetime of x- The former case is equivalent to having no additional terms due to interactions 
in ()28p and ()29p . and therefore after an initial build up, the condensate must fade quickly and will 
have no or very small contribution in the dark energy. In the latter case the contribution can be yet 
significant. 



In the approximation explained here, fc-dependent part of the integral (j56p is factorized and is pro- 
portional to F{k) defined in (j57p . It depends on the dimensionless variable: 

The contribution of fields having \[5zq\ < 1 grows exponentially as K becomes larger- most probably 
due to short range quantum fluctuations of the expectation values, but this behaviour is transient and 
eventually arrives to saturation. On the other hand, for fields with I/JzqI 2> 1, F{K) cx {c+ K"^)"^ . For 
large K (small distances) F[h) decreases with an slope index of —2, and for small K (large distances) 
it approaches a constant. At late times all the fields are in this regime and therefore this presents 
the late time behaviour of the expectation values. Therefore fc-independence is consistent with the 
observations of the dark energy at large scales. 

Finally, we note that the time dependence of (156p is determined by m/ Hq. This is similar to decoher- 
ence condition for massive scalar fields during inflation [25]. Only fields with small m/Ho can decohere 
and make classical density perturbations. This is another evidence for close relation between these 
processes as suggested in Appendix lAl 

3.2 Evolution of the condensate 

Now we use the result of the previous sections to estimate the evolution of the classical field x with 
cosmic time. Neglecting the self-interaction term in ()28p and (|29p . the general solutions of these 
equations at the WKB approximation level are: 

rv 

Xiv) = XiDq{a<i>r]) + X2Dq{-ais,ri) + g dr]' {Xk'^)G^{ri,ri') For decay mode (a) (59) 
2g dr]'{TXk)G^{r],r]') For decay mode (6) (60) 



where xij X2, Xi and X2 integration constants. We assume the following physically motivated 
initial conditions for these solutions: 

xivo) = x'ivo) = (61) 



The associated homogeneous equation of the differential equation (j28p is the same as equation (|34 
Therefore the Green's function G^{r],r]') is the same as the Feynman propagator (j39p . In the case 
of mode (b) we have added the WKB approximation to the homogeneous solutions. This correction 
is first order in the coupling g, thus we can use the Green's function ([39|) for this mode too. A 



more precise solution can be obtained by replacing Dg{a^7]) terms in ([39]) with the WKB corrected 
homogeneous solution Dg{—a^ J dr]{l — ^^^)*)- 

Using the asymptotic behaviour of Dq{z) at large z, we find that for both decay modes (a) and (b) 
the homogeneous part of the solutions ()59p and (j60|) is proportional to rj^^/"^ and therefore deceases 
with time. In the same way we can find time dependence of the special solution. For large rj it is 
proportional to r/^"*"^ for both decay modes. Here e is added by hand to parametrize the unknown time 
variation of the cosmological state |^'). If the mass of ^ is very small such that ~ 0, the solution 
of the evolution equations is approximately a plane wave in rj coordinate and x rf'^'' for the mode 
(b). For mode (a) this approximation does not modify the asymptotic behaviour of x- 

Finally the classical field 93 in comoving coordinates is x = o,^ ^ind varies as: 

if oc i~i/2+2+<: Assuming m 7^ for all fields, both decay modes. (62) 
If oc t-V2+3+e m$ « and mode (b). (63) 

For all cases, the late time behaviour of ip depends on e i.e. the back reaction. Based on ()41l) and 
qualitative arguments, we expect that e depends on the density and lifetime of X field, and therefore 
on the expansion of the Universe. If we want that at late times ip ^ cte., the only consistent model for 
the quintessence field is a model in which m$ 7^ but small. In such a model the condensate grows 
during the radiation dominated era but stops growing when matter becomes dominant. However, 
self-interaction should somehow modify this result. None the less, the strong limit on the clustering 
of the dark energy shows that self-interaction of the quintessence field can not be very strong and 
therefore its over all effect must be small. This simple model considered here favors particle physics 
models with a PNGB scalar field and a cyclical potential to protect the small mass of $ |24j . 

At late times when the energy density of the dark matter X becomes comparable to the density of the 
dark energy ~ l/2m|,(/?^, there is a strong feedback between expansion rate and the density of dark 
matter. Faster expansion rates bring down the density of the dark matter and therefore the rate of 93 
production decreases. This reduces the density of the dark energy and the expansion. We have seen 
the same feedback in classical treatment of this process in Ref. [5]. 

4 Outline 

Although we have not yet observed any elementary scalar field, we believe they play important roles in 
the foundation of fundamental forces in the nature, in Standard Model and in all suggested extension 
of it. The detailed studies of their behaviour and their consequences for other phenomena are however 
hampered by the fact that they can have complex non-linear self-interaction and interaction with other 
fields which make analytical calculations very difficult. 

We used quantum field theory techniques to determine the evolution of the classical component - the 
condensate - of a scalar field produced during the decay of a much heavier particle. Such a process 
had necessarily happened during the reheating of the Universe, and similar phenomena can happen 
- both as decay and as interaction - after the reheating if a scalar has interaction with other fields. 
We showed that a significant amount of condensate can only be obtained for light fields with masses 
comparable or smaller than Hubble constant at the initial time. By considering two decay modes, we 
also showed that the type of interactions has very important role in the cosmological evolution of the 
condensate and its contribution to dark energy. Therefore, the details of such a model depends on 
the particle physics of the fields. The general behaviour is nonetheless universal. Due to the coupling 
between quantum phenomena, i.e the micro-physics of the Universe, and its classical macroscopic 
content, obtaining an exact analytical solution is impossible. We left this task for future and by 
using simplifying assumptions we deduced the general aspects of the asymptotic time evolution of the 
condensate. We showed that one of the most important ingredient of this model is the unknown back- 
reaction of the decay and other physical processes on the Fock space of the Universe. We parametrized 



the time dependence of the Fock space with just one exponent oc 77^. If during matter dominated epoch 
~^ ^ £ ^ the condensate density evolves very slowly with time and has a behaviour similar to 
the observed dark energy. This range of epsilon is for minimal model without self-interaction for 
the scalar field. In a realistic context, both self-interaction and interaction with other field, specially 
Standard Model fields must be considered and they can change this range. 

The main goal of this work has been studying the contribution of the classical component of a gradually 
built-up field in the dark energy, specially in the context of a metastable heavy particle. Nonetheless, 
most of the results obtained here are applicable to other contexts and other models. Specifically, 
extensions of the Standard Model such as supersymmetric and supergravity models and string theory 
contain a large number of scalar fields. Condensation of these fields and their evolution and thereby 
their parameter space and symmetries can be constrained by the observation of the equation of state 
of the Universe. As we have shown here, the condensate evolution depends on a number of parameters 
including mass, decay rate of the producing particle, interactions and couplings. This is an additional 
and more detailed information with respect to simple density constraint for the scalar fields of particle 
physics models. 



Appendixes 

A A note on the relation between decoherence and condensation 

Using canonical representation and Bogolubov transformation for a free scalar field, one can find the 
relation between creation and annihilation operators at two different cosmic times: 

Wkiv)) Wkiv,Vo) vl{rj,rjo)J V4(%)/ 

^(??) = , '° ^=(4(r?) +afc(7?)) 

where a|.(%) and afc(%) are respectively creation and annihilation operator at an initial conformal time 
r/o, and Ufc(r/, ffc(??, and their conjugates are proportional to the solutions of the field equation at 
rj and 770 [22] • Using the definition of number operator N = a\.ak, it is easy to show that even without 
any self-interaction the number of particles at rj and r/o are not equal and cosmic expansion leads to 
particle production. The decoherence of these particles however needs an interaction to couple modes 
or fields [10] [29]. Although in a curved spacetime, specially during inflation, squeezed states can be 
achieved [27], the decoherence is not complete unless an interaction breaks the entanglement between 
degenerate quantum states by coupling causally unrelated modes - modes inside particle horizon to 
ones outside [28] [29] . 

We observe that for a free field in which modes are independent, the linearity of (164p and the properties 
of the creation/annihilation operators leads to (0|$(7/)|0) = 0, where |0) is the vacuum at r/o, i.e. there 
is no condensation. Therefore decoherence is a necessary but not sufficient condition for the formation 
of a condensate. 



(64) 
(65) 



B Free field Green's function on non- vacuum states 



In canonical representation, a free scalar field (j) can be decomposed to creation and annihilation 
operators of an orthogonal basis of the Fock space: 



(x) = ^Uk{x)ak+Ul{x)a\ , [ak,a\,] = 5kk' [afc,afc'] = [4,4,] = 



(66) 



where Uk{x) = Uk{ri)e ''^^ is a solution of the free field equation (f34|) . Quantization of 4> imposes the 
following relation: 

Uk{v,x)hl'C{r],y)-hl',{r],x)Ul{v,y) = i5^^Hx - y) (67) 

A Fock state 1^*) is constructed by multiple applications of the creation operator a\ on the vacuum 
state |0) defined by: 



afc|0)=0, V/c , |A:iA;2 . . . A;„) = 4i42 ■ • •4„|0) 

klk2...kn 

Applying these decompositions to 2-point free Green's function of (p, it can be written as: 

iGF{x,y) = {^\T(l){x)(l){y)\^) = 
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(70) 



From ()70p we can extract the expression for future and past propagators: 

iG>{x,y) ^ {^<p{x)<p{y)\^) = Y,Y. E ^kkS^k,k2...kfK{x)Uu{y) + 

k i k\k2...kn 

E 1 + E E 5kkS'^k,k2...kf Uk{x)Ul{y) 

k i kik2...kn 

iG<{x,y) ^ (^|0(x)0(y)|^) = ^^ ^ 5kk^^k,k2...kj'mxp*kiy) + 
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(72) 



Therefore, for free fields, Gpix^y) on any density operator p can be written as a linear expansion 
with respect to Uk{x)U^{y) and IJl{x)Uk{y), the independent solutions of the free field equation. The 
contribution from a non- vacuum state |^) appears in the coefficients of the expansion, just as the ini- 
tial/boundary condition effects appear in the expansion coefficients of the propagator (j39p . Projection 
coefficients \'^kik2...k„\'^ determine the momentum distribution of the states in the environment. In 
the simplest case where quantum correlation between particles is negligible, it is proportional to one 
particle momentum distribution f{k) that in the case of a thermal environment is Boltzmann, Fermi 
or Bose-Einstein distribution. 

Here a general comment about expectation values including 2-point Green's function is in order. An 
expectation value i.e. an inner product in the Fock space of a quantum system, this includes operators 
at different spacetime points is meaningful only if these points share the same or isomorphic Fock 
spaces. The state |\I') in (j70|) for which the expectation value is calculated must be a member of 
the Fock space at x and y. However, in the fast changing early Universe, it is not evident that such 



a condition exists unless x and y are enough close to each other. This reflects the back-reaction 
of interactions during and after reheating on the Fock space and coupling between operators and 
states as mentioned in Sec l2.31 Therefore, |^') has an implicit spacetime dependence and one-particle 
distribution function gets the familiar form f(x, k) used in the construction of classical Boltzmann 
equation [26] . 



C Propagators and evolution in matter dominated epoch 

In the matter dominated epoch the relation between comoving and conformal time is defined as: 

fdt 3to 
??=/— = %(-) , Vo = — (73) 
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By applying (j74p to the field equations (|28p and (j29p and neglecting interactions, the field equation 
for both modes is the following: 

x" + {k' + ^^^--,)x = (75) 
Vo T 

where x presents one of the fields X) X or A. For two special cases of m = and A;^ = this equation 
has exact analytical solutions: 

1 VvJ±l{k''l) Fot m = 

As usual, in each case the WKB approximation can be used to find an approximate solution when the 
neglected terms are not zero. Another possible approximation is a linear interpolation between two 
cases in (1761). 



Bessel functions J±i/2 ^-iid J±z/2 have close analytical expressions: 



and 



xiv) 



Ji{x) = \ — sinx , J_i{x) = \ — cos a; (77) 
2 V vrx 2 V vrx 

J-i{x) = -Ji{x) - J^i{x) , J_-i{x) = -Ji{x) --J^i{x) (78) 

2X2 2 2 2X2 

Xi sin{j3rf) + X2 cos(/37y^)^ For /c^ = 



- cos(A:r/)) + ^^{- sin(fcr?) - ^^)) For m = 



(79) 



where xi) X2) x'l &iid X2 integration constants. At late times and for large scales (small k) the 
mass term in (|75p is dominant and /c^ = approximation can be applied. Using (j79p . we find that 
integrals analog to (|52p (replacing with x from (j76p ) have a late time behaviour oc rf^"^^ where 
N is the number of fields in the expectation value. If we assume that the mass of <I> is very small and 
m = 0, this approximation can be applied to this field. Under these conditions the expectation values 
containing one x field are oc rf""^^^"^) where here N is the number of other fields. Finally from these 
results and the general solutions (|59p and (|60p (again replacing lA}^ with x) we conclude that: 

Lp OC , For m 7^ 0, both modes (80) 

4 + 1! 

if (xt 3 , For m$ = and mode (b) (81) 



Comparing (j80|) and (jSip with (j62p . and (j63p and assuming the same e, we find that in the matter 
dominated epoch the time evolution is slower than radiation dominated epoch. In the matter dom- 
inated epoch the expansion of the Universe is faster, the decay of X particles is slower because for 
the same average density the Universe is younger, and therefore less X particles have decayed. This 
compensates the density reduction due to a faster expansion. Similar to the radiation dominated 
epoch, the final time evolution rate depends on the unknown parameter e. 
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